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I. INTRODUCTION 



Numerical studies, from variational calculation to exact diagonalization, all indicate that the 

quasiparticle generated by introducing one hole into a two-dimensional quantum antiferromagnet 

has the same nature as a string state in the t — J z model. Based on this observation, we attempt 

to visualize the quasiparticle formation and subsequent coherent propagation at low energy by 

studying the generalized t — J z ~ J± model in which we first diagonalize the t — J z model and then 

perform a degenerate perturbation in J±. We construct the quasiparticle state and derive an effective 

Hamiltonian describing the coherent propagation of the quasiparticle and its interaction with the 

spin wave excitations in the presence of the Neel order. We expect that qualitative properties of 

the quasiparticle remain intact when analytically continuing J± from the anisotropic J± < J z to 

the isotropic J± — J z limit, despite the fact that the spin wave excitations change from gapful to 

gapless. Extrapolating to J± — J z , our quasiparticle dispersion and spectral weight compare well 

^ with the exact numerical results for small clusters. 
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Since the beginning of the era of high temperature superconductivity, it has been suggested that the one-band 
Hubbard model ffi^the closely related t — J model describing the CuC>2 plane may be the basic model for cuprate 
superconductors.Era While the early suggestion was based on the high energy properties of the cuprates, such as 
electronic structure and spectroscopy, it is remarkable that many peculiar-low energy properties of the cuprates have 
been reproduced in exact numerical studies of the Hubbard or t — J modclE Thus, understanding the two-dimensional 
Hubbard or t — J model stands out as a central problem in the field of high temperature superconductivity. 

Whileuthe Hubbard and t — J models at half- filling describing the insulating parent compound-are. well understood 
by now,Q present understanding of the effect of doping mostly derives from numerical studies .H-LnlJ Essentially, the 
numerical results can be summarized as follows: (i) Doping the two-dimensional (2D) antiferromagnet with one hole 
generates a quasiparticle with finite spectral weight which implies that the quasiparticle carries charge e and spin 
1/2. (ii) The quasiparticle dispersion has a bandwidth of order J, is almost flat near (it, 0) and (0, 7r), and the general 
structure agrees with the recent photoemission results, (iii) The remaining spectral weight is mostly distributed in 
k> ' the energy rangexprresponding to the mid-infrared absorption universally observed in the optical experiments of the 
high T c cuprates&j. (iv) There is tentative evideace-andicating that the quasiparticle band found in the one-hole 
case remains fairly robust under the finite doping.£3il3 The latter may have already received experimental support.tll 
In this paper, we attempt to provide an intuitive way to understand all numerical results, and derive a low-energy 
effective Hamiltonian for the one-hole problem. We shall also point out the extension to finite doping when the Neel 
order is destroyed. 

The exact numerical results on the one-hole problem in small clusters suggest an adiabatic picture for the doped 
antiferromagnet. For the t — J model with t > J, the hole hopping takes place at a typical time scale of 1/t while 
the time scale of spin fluctuations is at least 1/J("S> 1/t). When considering the hole hopping, it is reasonable to 
take a snapshot of the spin configuration and freeze it. When the system possesses either long range or short range 
antiferromagnetic correlation, a typical spin configuration also possesses the same kind of correlation. As the hole 
hops in two dimensions, it leaves behind a string of overturned spins along its trace. These overturned spins are 
energetically frustrated with their neighboring spins surrounding the trace. The energy increase due to the frustration 
is roughly proportional to the length of the string. Therefore, the hole is linearly confined to its original site. Due 
to this confinement, the hole is dressed by a cloud of local distortion of the spin configuration. The quasiparticle 
is just such a composite object. The formation of the quasiparticle does not require the existence of long range 
antiferromagnetic order. The quasiparticle should be stable as long as its size is smaller than the magnetic correlation 
length, ft is also important to realize that the formation of the quasiparticle is a high-energy process at the scale of the 
hopping amplitude t. Exact diagonalization study of small clusters provides accurate description of the quasiparticle 
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structure. At the lower-energy scale of order J, the physical properties of the doped antiferromagnet are described 
by the coherent propagation of the quasiparticle and its interaction with the spin wave excitations. It is at this 
lower-energy scale that an effective Hamiltonian is needed. 

Our approach amounts to a degenerate perturbative treatment of the transverse spin fluctuation term J± of the 
t— J z — Jj_ model. We first construct a basis of states by diagonalizing the t— J z model. Then we project the transverse 
J_l term onto the constructed basis of states and use the resuhmfcaiairix elements as an effective representation of 
the transverse spin fluctuation term. Previous analytical studiea3il3~E3 have concentrated on solving Bogoliubov- 
deGennes-type equations (summing up the noncrossing diagrams) in which the spin wave excitations are represented 
by Holstein-Primakoff bosonsu and the hole hopping term is treated as the interaction. Although the plausible 
argument for the "dominant-pole approximation" put forward by Kane, Lee, and Readlia indicated the existence of 
a low-energy quasiparticle, it has been found that a self-consistent solution is needed to obtain the correct featurpa 
of the quasiparticle dispersion and spectral weight. The self-consistent solution can be obtained only numericallyi3'c2l 
The advantage of our approach is the clear separation between the quasiparticle formation at high energy from the 
low-energy coherent propagation of the quasiparticle. The parameters in the effective Hamiltonian describing the 
quasiparticle propagation are evaluated directly, though approximately. Most importantly, the form of the effective 
Hamiltonian does not depend on the detailed calculations and is completely determined by the spin and charge 
quantum numbers of the quasiparticle. 

The specific nature of the quasiparticle, that it is essentially same as a string state as explicitly shown by Eder 
and Ohta,Ej implies that the qualitative properties of the quasiparticle remain intact when analytically continuing J± 
from the anisotropic J± < J z to the isotropic limit J± = J z , despite the fact that the spin wave excitations change 
from gapful to gapless. We have evaluated the parameters of the effective Hamiltonian to first order in J± which 
is appropriate when Jy_ <C J z . Nevertheless, by taking J± = J z our quasiparticle dispersion and spectral weight 
compare qualitatively and quantitatively well with the exact diagonalization results for the isotropic t — J model. The 
quantitative agreement is expected to improve should our results be compared with the exact numerical results for 
the anisotropic t — J z — Jj_ model if available. These expectations can be critically tested. 

Although the hole hopping amplitude t is much larger than J, the bare hole hopping is localized inside the quasi- 
particle and the associated excitations are gapful. The true low-energy charge excitations are those corresponding to 
the coherent propagation of the quasiparticle. The velocity of the quasiparticle near the bottom of the quasiparti- 
cle band is much smaller than the spin wave velocity. Therefore, the spin wave excitation is not dissolved into the 
quasiparticle-hole continuum excitations. As a result, the low-energy spin wave excitations cannot decay by exciting 
the quasinaxticle into a higher-energy state of the quasiparticle band. The linear dispersion of the spin excitations is 
preserved JEjlJ contrary to naive expectation based on the fact t~^> J. 

The layout of the paper is as follows. In Secjll], we derive the effective Hamiltonian for one hole in the antiferromag- 
net. In Sec. [II, we compare our quasiparticle dispersion and spectral weight with the exact diagonalization results in 
small clusters. We conclude the paper in Sec. IV. 



II. ADIABATIC DESCRIPTION OF A DOPED ANTIFERROMAGNET 



In this section, we shall derive an effective Hamiltonian describing the quasiparticle motion in a doped antiferro- 
magnet and the interaction of the quasiparticle with long wavelength spin wave excitations in the presence of Neel 
order. We shall first construct all low-energy eigenstates of the t — J z Hamiltonian. Then we represent the transverse 
Heisenberg interaction J± term by its matrix elements within the basis of these eigenstates of the t — J z Hamiltonian. 



A. Construction of the basis states 



The t — J model Hamiltonian has the form 



H = H + H 1) (1) 
H = -tJ2 [(1 - «,,-akW(l - ni,-«r) + h.c] + S t S h ( 2 ) 



<i.j> <ij> 



H t = J ± J2 ( S i S ? + S ? S j)> (3) 



<hj> 



where c\ a creates an electron at lattice site i with spin a, rii_ a — c^ ia Ci a is the electron density operator, and Si — 
c \n®nvCi V /2 is the spin operator. The summation <i,j> is limited to the nearest neighbor links. Unlike the Neel 
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state in the undoped case, the ground state of the t — J z Hamiltonian Hq in the presence of one hole has macroscopic 
near degeneracy. These ground states can be generated by removing a spin at an arbitrary site in the Neel state and 
allowing the hole to hop. As the hole hops, it leaves behind a string of overturned spins. Therefore, the hole suffers 
a linear binding potential to its original site. The eigenstates of the hole consist of discrete bound states with energy 
spacing uniquely characterized by an energy scale (J z /t) 2 / 3 t. For the purpose of deriving a low-energy effective theory, 
we only need to retain the ground state. . ., . 

There is an excellent approximation to obtain the ground states of the t — J z Hamiltonian.E3t3 The approximation 
amounts to neglecting winding paths (Bethe lattice approximation) . We briefly describe this approximation here with 
some minor improvement. We shall also use slightly more concrete notations to meet our needs. The perfect classical 
Neel state with only one hole can be labeled by the position of the hole, r. There are no frustrated spins in this state. 
Applying the hole hopping operator in H Q to this state generates four new configurations, each with a string of a 
length of one lattice spacing, |r, {a}), where a can be one of the four unit lattice vectors ±x or ±y, representing the 
hopping of the hole from r to r + a. From the four configurations we construct an "s-wave" state with string length 
1, 

|r, 1> = | J2 l r »>- ( 4 ) 

a— ±x,±y 

Other states with longer strings are constructed similarly, 

M)= 2x3 i_ 1)/2 Yl |r,{ai,a 2 ,--- a/ }), (5) 

ai—ai 

where the summation is implicitly subject to the constraint a. L ^ — a.;_i to prevent hole retracing, here and throughout 
the paper. The Bethe lattice approximation is to use the matrix elements within the set of states |r, I) to represent 
the t — J z Hamiltonian for the corresponding subspace, 

ff o |r,0) = J z |r,0)-2i|r,l), 

H \r, 1) - p z \r, 1) - It |r, 0) - VSt |r, 2), (6) 
H \t, I) ~ e/|r, I) - V3t (|r, I - 1) + |r, I + 1)) , I > 2, 

e l = J z {l,v\ Us; + ^\r,l). (7) 

<i,j> ^ ' 

The usage <Athe average energy e; slightly improves the approximation (||) for the t — J z Hamiltonian used by Shraiman 
and Siggia.EZI In this paper, we shall directly evaluate the first few values of e; up to I = 7. For larger I values, we shall 
use the extrapolation of e; from small /. The values of e; used in this paper are listed in Table Q. The extrapolation 
is e; = 2.574 — 0.0156(— 1)' — 0.7857// + 0.6857L Numerically diagonalizing (||) determines the coefficients ui in the 
approximate ground state of the t — J z Hamiltonian, 

oo 

H \r) = -E \r), |r) = |r,l). (8) 

1=0 

There is a macroscopic number of such ground states since we can choose r to be any lattice site. The size of |r) 
depends on J z /t ratio: The smaller the ratio, the bigger the size. 

In the two-dimensional square lattice, the macroscopic degeneracy of |r) is removed by the winding paths (so-called 
Trugman process) .0 This is reflected in the fact that the approximate ground states |r) labeled by r and given by 
Eq. (||) are not strictly orthogonal to each other. In the presence of Neel order, we can divide the two-dimensional 
square lattice into A and B sublattices. For r belonging to the same sublattice, there is a small overlap. An example 
is shown in Fig. [y. This configuration belongs to both |r = ro, I = 5) and |r = r2, 1 = 1) according to the prescription 
(||), causing an overlap between |ro) and |r2) states. For J z /t = 0.3, the overlap between |r) and |r + x + y) is 0.0409. 
The overlap between |r) and |r + 2x) is 0.0019, negligible even compared to (r|r + x + y). For r belonging to different 
sublattices, any two states are orthogonal since they have different S^ ot = J^i $1 ■ The small corrections due to the 
nonorthogonality of |r) will be taken into account later by using tight-binding approximation to study the propagation 
of |r) as a whole. 



where 
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The set of states |r) is far from complete. One way to complement the set of states |r) is to include "non-s-wave" 
string states. Instead ([|), for instance, we can form the d state, 

|r, 1)„ = \ (|r, {x}) - |r, {y}) + |r, {-x}) - |r, {-y})) . (9) 

Applying the hole hopping operator in Hq to (|^) will generate a new set of states. Diagonalizing the t — J z Hamiltonian 
within the corresponding new subspace will give us new eigenstates of the t — J z Hamiltonian whose energies are about 
1.5 J z higher than the states |r) given by Eq. (||). Obviously, the new eigenstates are also localized around r. For the 
same reason — the existence of the winding paths — the new eigenstates centered at r are not strictly orthogonal to 
the "s-wave" states |r') centered at a different site r' in the same sublattice and not too far away from r. We note 
that if r = r' there is no overlap due to the different rotational symmetries. The spin fluctuation term (|J) also has 
nonvanishing matrix elements between the "s-wave" states |r') and the new eigenstates. However, both overlaps and 
matrix elements are very small. Thus, these new eigenstates of the t — J z Hamiltonian generated from "non-s-wave" 
string configurations are almost disconnected from the subspace where the states |r) "live". The "non-s-wave" states 
will be omitted in our low-energy effective theory. 

The eigenstates of the t— J z Hamiltonian that we have discussed so far contain overturned spins only along the trace 
of the hole. If we neglect the winding paths, the overturned spins form a continuous string in each spin configuration. 
We need to complement the set of states |r) by including one or more isolated overturned spins outside the string. 
Inside the eigenstates |r), the hole is localized around the site r, and so are all the strings of overturned spins. Thus, 
if the isolated overturned spins are far from r, they will not interfere with the strings. In this case, we have simple 
superposition of the isolated overturned spins and the string state |r). More delicate situations occur when the isolated 
overturned spins are located close to the center site r of the string state. Our goal is to find all new eigenstates of the 
t — J z Hamiltonian that can be connected to the string states |r) by applying Hi. 

We start from a configuration with one isolated overturned spin at r s and a hole at r in an otherwise perfect 
Neel background. We demand that a new subspace be generated by applying the hole hopping operator to this 
configuration and new eigenstates of the t — J z Hamiltonian be generated by diagonalizing the t — J z Hamiltonian 
within this new subspace. In order for H\ to have nonvanishing matrix elements between the new states and the 
old states given by Eq. (JsJ) , the site of the hole r and the site of the overturned spin r s must belong to different 
sublattices. Otherwise, the generated new state has S* ot = ±3/2 and cannot be connected to the states (||) by Hi 
since the states (||) have Sf ot = ±1/2. Since we require r s to be close to r, the first possibility is for r s to be a nearest 
neighbor of r. But this is a string of length 1 |r s , {r — r s }) which we have already considered. Thus, no new subspace 
will be generated from this configuration. The next possibility is for r s to be farther away from r: r s — r = 2x + y 
or r s — r = 3x. Letting the hole hop in these configurations will generate new subspaces. In principle, we have to 
diagonalize the t — J z Hamiltonian within these new subspaces and keep the low-energy eigenstates. In practice, we 
shall do it approximately as in (||) . In complete analogy with (|^) , we construct generalized states 

Vs,r,l) = - - 2(1-1) /2 l r ^r,{ai,a 2 ,---a ( }). (10) 

ai-'-a; 

Then we use the matrix elements of the t — J z Hamiltonian within the set of states |r s , r, I) for I = 0, 1, • • • oo as an 
approximate representation of the t — J z Hamiltonian. This representation has exactly the same form as ([|) except 
for different diagonal energies e[: 

ff |r s ,r, 0) = 3J z |r s , r, 0) - 2i |r„ r, 1), 

H \v a ,v,l) = ^J z |r s ,r,l) -2t|r s ,r,0> - V3t|r s ,r,2), (11) 
flo|r„r,J> ~e;|r s ,r,0-V3i (|r s ,r, I - 1) + \r s ,r,l + 1)) , I > 2, 

where 

e[ = J z {l,v,v s \ (s?S* + l)\r s ,r,l). (12) 

<i,j> > ' 

The values of ej used in this paper are also listed in Table fl. For large I, we use extrapolated values from ej listed in 
Table | in the same way as we did for e; . Diagonalizing ( |Tl| ) gives us new approximate eigenstates 

H \r s ,r) = -E' Q (r s -r)|r s ,r), 

OO 

|r. l r)=£u{(r,-r)|r a> r J l). (13) 

1=0 
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Following the same line of reasoning, we can construct further approximate eigenstates of the t — J z Hamiltonian. 
A straightforward extension is to consider r s being further away from r, for instance, r s — r = 5x, 4x + y, or 3x + 2y. 
Since in the realistic situation J/t — 0.2 ~ 0.4, the size of the string state |r) is limited within I < 5 inside (||). 
Therefore, we can approximately treat the isolated overturned spin at r s as independent of r if their separation is 
five or more lattice spacings. We should also try to generate new eigenstates by applying the hole hopping operator 
in Hq to configurations with one hole and an isolated pair of neighboring overturned spins in an otherwise perfect 
Neel background. If the hole is a nearest neighbor of the pair of overturned spins, this is a string of length 2 that 
we have already considered in (|^) . As the next possibility, the hole could be separated from the pair of overturned 
spins by one lattice site. One such example is shown in the upper panel of Fig. |^. However, this configuration again 
has already appeared in ( |l(i| ) in the I = 3 state as illustrated in the lower panel of Fig. |[ No new eigenstates of the 
t — J z Hamiltonian can be generated from this configuration. Thus, we only need to consider configurations in which 
the hole is separated from the pair of overturned spins by two or more lattice sites. In this paper we shall neglect the 
correlation between the hole hopping and the pair of overturned spins in this kind of configurations. Later we shall 
see that what we have neglected are the interaction vertices with four operators, two quasiparticle and two magnon 
operators. 



B. Derivation of the effective Hamiltonian 



To obtain the low-energy effective Hamiltonian, we identify the many-body state |r) with a quasiparticle located at 
site r. It is obvious that the quasiparticle carries charge e and spin 1/2 as a hole because there is local deficiency of 
both charge and spin. Furthermore, the quasiparticle spin points in the opposite directions for r belonging to different 
sublattices. Thus we introduce a fermionic operator fl a to create the state |r) from the Neel background. In the Neel 
state, we divide the square lattice into A and B sublattices and assume that the spins in the A sublattice point in the 
up direction. If r € A, then the quasiparticle spin is down. We shall approximate |r s ,r) ~ 5+/^ J Neel) for r G A 

and r s G B. Using the Holstein-Primakoff representation of the spin operator ,0 — b\ 1 for i G A and Sf — b\ 2 for 
i £ 6, we can write 

|r s , r) ~ bl jljmel) for r G A, r s G B. (14) 

The expression for r G B and r s G A is similar. The correlation between the hole hopping and the overturned spin at 
r s lowers the energy of the state |r s ,r) with respect to that of two independent &J s and /J . This energy difference 
amounts to an interaction vertex of the type /J a f r ,ab\ b r ^. In this paper, we shall neglect all interaction vertices with 
four operators and assume that their effects can be properly accounted for by small renormalization of the parameters 
in the effective Hamiltonian. This is supported by the fact that quantum fluctuations only generate a small density 
of overturned spins in the ground state of the 2D Heisenberg antiferromagnet. 

We have pointed out that the essential difference from the undoped case lies in the macroscopic near degeneracy of 
the ground states of the t—J z Hamiltonian. The transverse spin fluctuation J_l term in the t — J z — J± Hamiltonian 
(§) lifts this dc generacy and generates quasiparticle dispersion. Therefore, only degenerate perturbation in J± is 
permitted even for J± <C J z . In the following we shall represent the transverse term (g) by its matrix elements within 
the basis of the eigenstates of the t — J z Hamiltonian that we have constructed in the previous subsection. A typical 
example illustrating how the spin fluctuations generate quasiparticle hopping is shown in Fig. |^. Exchanging the first 
two spins on a string cuts off its length by two lattice units. Since the location of the quasiparticle is determined by 
the common starting point of the strings inside @ and (^) , this process moves the quasiparticle to a nearby site in 
the same sublattice. Similarly, exchanging two spins near the starting point but outside a string increases the string 
length by two lattice units and also generates quasiparticle hopping. We note that a quasiparticle cannot hop to the 
other sublattice simply because the quasiparticle spin directions are opposite in different sublattices. To obtain the 
effective hopping amplitude of the quasiparticle, we construct the propagating quasiparticle state 

|k) = ^=£exp(ik.r)|r), (15) 

where 1 / y/Nk is-the normalization factor. Then we calculate the quasiparticle dispersion according to the tight-binding 
approximation^ 

e k =^<k|ff|k>, (16) 
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where H is the full t — J Hamiltonian. Using the fact that the string state |r) is an exact eigenstate of the Hamiltonian 



H r = -t 2^ [(l-n it - a )cl a c ja {l-n jt - a )+ h.c] + 2^e/|r,Z)(r,Z|, (17) 

<i,j> I 

we can rewrite 

e k = const + — exp(ikR)(r + R\H - H r \r). (18) 
k ReA 

The normalization factor can be expanded in a power series, 

— — ~ 1 — 4(r + x + y|r) cos k x cosk y . (19) 

For J z /t = 0.2, 0.3, 0.4, 0.5, and 0.6, we have found (r + x + y|r) = 0.065, 0.041, 0.027, 0.018, and 0.012 respectively. 
The direct overlap (r + 2x|r) is neglected in Jie| ) because it is too small even compared to (r + x + y|r). Similarly, 

the matrix element (r + 2x| ^ J z j> S*Sj — J2i £ il r > 0( r : ^l) l r )> which appears in the calculation of Eq. (|is|), can 
also be neglected because the nonvanishing contributions to this matrix element have the same origin as in the direct 
overlap (r + 2x|r). After some straightforward algebra, we obtain from Eq. ( |i~8| ) 

ek = const + Aa% cos k x cos k y + 2a2(cos 2k x + cos 2k y ) + C(cos k x cos 3k y , cos 2k x cos 2k y ), (20) 

where the high order trigonometric functions corresponding to longer range quasiparticlc hopping arc neglected. The 
coefficients in the dispersion are given by 

ai = (r + x + y\H - H r \r) - [(r|ffi|r) + 2(r + 2x|# 1 |r)] (r + x + y|r), (21) 
a 2 = (r + 2x\H 1 \r)-2{r + x + y\r){r + x + y\H-H r \r). (22) 

This quasiparticle dispersion can be described by a tight-binding Hamiltonian 



*tA E 



2 

reA,u=± 



ai 22 + 012 

a=±x±y a=±2x,±2y. 



(fljr+*,o + h.c.) , (23) 



where r G A denotes summation over only one sublattice. The summation over A and B sublattices is converted into 
the spin index a summation. 

To obtain the coefficients a\ and we need to calculate the matrix elements of the type (r2|0|ri) for several 
different operators O. We employ the following approximation when calculating these matrix elements: (r2|0|ri) = 
J2i 1 +i 2 <io u h u h ( r 2 , ^2 |0|ri , Zi) . Each matrix element in the sum is then computed straightforwardly. By investigating 
how fast the eigenvector ui decreases as a function of the string length, one can be convinced that the approximation 
^i + h < 10 is sufficient for J z /t > 0.2. The values of a\ and 02 for different J z /t ratios are listed in Table [Dj. 

The spin fluctuation term Hi exchanges two neighboring spins. For most cases, this pair of exchanged spins is far 
away from the hole. If so, the situation is completely same as in the undoped case. Thus, the Hamiltonian describing 
long wavelength spin excitations using the Holstein-Primakoff bosons has the familiar form,Eij upon neglecting the 
spin wave interactions, 



Hh — 2J Z 



E (t>tA + *,2 + Kib W ). (24) 



2 

ieA a=±x,±y 



In principle, the summation over the nearest neighbors in (|24| ) should avoid a small region surrounding the hole to 
prevent double counting. However, this restriction only results in interaction vertices involving four operators of the 
type P fb^b or f'fb'b'. These interactions are neglected in this paper. Physically, these interactions occur because 
the quasiparticle has a finite spatial extent of a few lattice spacing. But for long wavelength spin excitations the 
quasiparticle can be viewed as a point particle. 

The spin fluctuation term Hi also connects the states given by Eqs. (^) and (|l^). A typical example is shown in 
Fig. § This process allows the quasiparticlc to jump over two sites, leaving behind an isolated overturned spin. Wc 
directly evaluate the following two overlap coefficients: 
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Ai(J 2 /t) = — (r|ff 1 |r,r + 2x + y), 

X 2 (J g /t) = ^-(T\H 1 \T,T + 3x). 



(25) 
(26) 



The nonorthogonality of the approximate eigenstates |r) has a negligible effect on the coefficients Ai and A2. The 
results are also listed in Table ||. The corresponding part of the effective Hamiltonian is, using Eq. (|l4]), 



H f b - Jl - E 



Ai 



E +^ E 



r£A L a=±2x±y,±x±2y a=±3x,±3y. 



^E 



a, y: +>» e 

a=±2x±y,±x±2y a=±3x,±3y. 



(27) 



The complete effective Hamiltonian is, upon neglecting magnon interactions and quasiparticle-magnon interactions 
containing four or more operators, 



H e s = Hf + Hi, + H 



fb- 



(28) 
sum- 



Since inside H c g the lattice site summation only extends within one sublattice, the corresponding momentum . 
mation should be limited within the antiferromagnetic Brillouin zone(AFBZ). This is due to the presence of the Neel 
order which breaks the symmetry between the two sublattices. In the momentum space, the three parts of H c g have 
the following forms: 



H f = E £ k fLh, 

keAFBZ,(T=± 



H b = Y [ 2 -^( 6 q,l & q,l + & q,2 & q,2) + J± (COS q x + COSq y ) f& q)1 &l q)2 + frq4 fo -q,: 



qeAFBZ 



H 



fb 



4ff E * 



(29) 
(30) 
(31) 



/k, T /k+qa & ql + /k,l/k+q,T & q2 + h - C 

k, qeAFBZ 

where M is the number of sites of one sublattice. The quasiparticle-magnon interaction vertex is given by 

Ak = 2 J± [2Ai (cos 2k x cos k y + cos k x cos 2k y ) + A2 (cos 3k x + cos 3k y )] . (32) 
The dispersion ek is given by Eq. (p0|). 



III. QUASIPARTICLE PROPERTIES 

The quasiparticle dispersion parametrized by a\ and «2(see TaJaki,p)_£orrectly captures all qualitative features 
found in the numerical studies of the t — J model with one hole.0l3BBEl These include (i) the right magnitude of 
band width given by 4(«i + 2a2) (ii) the tight-binding dispersion of the same-sublattice hopping type with the bottom 
of the band at (%/2,tt/2), and (hi) the small energy difference between the quasiparticle energies at the momenta 
(tt/2, 7r/2) and (ir,0). This small energy difference results from the closeness between the values a% and 2o!2- In 
our approach, this is due to the fact that a quasiparticle can hop from the site r to r + x + y by exchanging spins 
either between the sites r and r + x or between r and r + y. In contrast, the two spins on the sites r and r + x 
must be exchanged in order to move the quasiparticle from r to r + 2x as illustrated in Fig. |^. Taking t = 1 and 
Jj_ = J z = 0.4 as an example, our dispersion is given by £k = 0.387 cos k x cosk y + 0.110(co&2At; + cos2/cj,). This is to 
be compared with ek = 0.34 cos k x cos k y + 0.13(cos2fc a; + cos 2k y ) reported in the literature.Hl23 From Table || and by 
setting J± — 0, we notice that the quasiparticle has a small dispersion even, for the t — J z model and the minimum of 
the band is located at k = (0, 0), in agreement with the previous result.EHlJ 

The quasiparticle spectral weight is determined by the overlap between a propagating bare hole and the quasiparticle 
wave function |k). From Eqs. ( |l5| ) and (0), it is simply given by Uq/N^. For different t/J z ratios, the quasiparticle 
spectral weight at the bottom of the band (7r/2, tt/2) is listed in Table Within our first order degenerate perturbation 
in Ji, the effect of the spin fluctuations on the internal structure of the quasiparticle is not included. The consequence 
is that the quasiparticle spectral weight only depends on J z /t. In the exact diagonalization study of the isotropic 
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t — J model, the spectral weight is momentum dependent. Generally, the spectral weight is found smaller for the 
high-energy states of the quasiparticle dispersion. This reshuffling of the spectral weight could be the reason for the 
systematic underestimation of the spectral weight in our first order degenerate perturbation. Physically, one can 
imagine that the quantum spin fluctuations smear the energy gap between the localized string states of the t — J z 
model and spread the spectral weight of the first excited state of the t — J z model toward the low energy side. If the 
nature of the quasiparticle is not changed when including the J± effect on the quasiparticle structure, one expects 
to improve the first order degenerate perturbation by optimizing the coefficients ui in Eq. (j^). Eder and Ohta have 
achieved this by constructing a quasiparticle operate* in which the coefficients corresponding to u; are optimized 
against their exact diagonalization spectral function.cj Nevertheless, our first order degenerate perturbation results 
give a fairly good estimate of the quasiparticle spectral weight at the ground state wave vector (tt/2,it/2) of the 
isotropic t — J modeLtrB'Eij 

Besides the effect of the quantum spin fluctuations on the quasiparticle structure, the quasiparticle motion also 
suffers scattering from the spin wave excitations due to the interaction terms in -H/b. However, the second order 
corrections to the quasiparticle dispersion and spectral weight usually does not exceed a few percent if we use the 
values Ai and A2 listed in Table ||. This is most due to the fact that the vertex Ak given by ( |32|) vanishes around the 
bottom of the band. To the accuracy of applying the first order degenerate perturbation to the isotropic t — J model, 
these corrections can be neglected. Strictly speaking, the spin wave excitations can mediate long range interactions 
between quasiparticles in the presence of Neel order. However, for the interaction vertex given by Ak the disorder in 
real materials should localize the quasiparticles before the effect of the long range interactions becomes significant. 
For completeness, we include the expression for the second order self-energy of the quasiparticle here. For the one-hole 
problem with J z = Jx = J ', we found 

. , 1 + JA k+q [2A k + q - A k (cosg x + cosg y )]/w q 

Z(k,lUJ n ) = — -. > ^ ; , (33) 

where u> n = (2n + \)nT is the Matsubara frequency and w q = 2 J \J\ — (cos q x 4- cos q y ) 2 /4. We note that for the 
one-hole problem the chemical potential n is chosen such that ek > for all k. 



IV. CONCLUSION 



We have developed a degenerate perturbative treatment of the t — J z — Jj_ model based on the observation that 
the quasiparticle has the. same nature as a string state. Evidence supporiiiLg this picture is abundant from early 
variational calculationsOo to the more recent exact diagonalization studiesoEl A prerequisite for the success of the 
variational calculation is that the elementary excitations can be described by a hole dressed with strings of overturned 
spins. Recently, Eder and Ohta calculated the one-particle spectral function .using a composite operator precisely 
describing a string state |r) given by (|8|) except with optimized coefficients it;.E3 They found that the quasiparticle 
peak is greatly enhanced while the incoherent part is suppressed when the composite operator is used instead of 
the bare hole operator. This unambiguously reveals the nature of the quasiparticle since by progressively improving 
the composite operator one will presumably reach a point where only the quasiparticle peak is left. Thus, an exact 
quasiparticle operator can be constructed at least conceptually. I— ■ 

Is the quasiparticle discussed in this paper related to the so-called spin bag_ introduced by Schrieffer et alEl? We 
will postpone the detailed discussion of this question to another publication.^ The original spin bag was introduced 
in connection with the Hubbard model, and the bag is a local reduction of the Neel order parameter associated with 
the bare hole. The size of the spin bag is reduced when U is increased, whereas the quasiparticle of the t — J model, 
discussed here, is increasing in size, when J is decreased, i.e., when U is increasing. In fact, we shall argue, that in 
the limit of large U Hubbard model the spin bag will become thej-hare hole in the t — J model. We shall present a 
unified treatment of the Hubbard model in the future publication.El 

With finite doping, it is known that the long range antiferromagnetic order is destroyed. If the quasiparticle is still 
formed and stable, the propagation of the quasiparticle will suffer frustration at the length scale of the magnetic cor- 
relation length because of the fluctuations of the local spin directions. Traveling along a closed path, the quasiparticle 
will pick up a Berry phase, resulting in a loss of phase coherence. These issues will also be discussed in the future 
publication. 
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TABLE I. Average energies of the string configurations of different lengths 



1 




e[/ J z for r s — r = 2x + y 


eJ/J 2 for r s — r = 3x 





1 


3 


3 


1 


2.5 


4.5 


4.5 


2 


3.5 


5.375 


5.4583 


3 


4.3889 


6.0139 


6.2639 


4 


5.1296 


6.7454 


6.963 


5 


5.8704 


7.412 


7.6142 


6 


6.5535 


8.0828 


8.286 


7 


7.2531 







TABLE II. The parameters of effective Hamiltonian evaluated using the string state approximation. Zh(k) is the quasi- 
particle spectral weight. The values of Zh in the parentheses-are the exact diagoimlization results of the 26 site cluster at the 
ground state wave vector k* for the t — J model from RefsBI— I. Note that, in Refi— 1, Zh increases with the cluster size for the 
t — J z model at J z /t = 0.3 while Zh at k* decreases with the cluster size for the isotropic t ~ J model at J/t = 0.3. 



Jz/t Z h (tt/2, ty/2) oi «2 Al A2_ 

0.2 0.160(0.213) 0.305 J ± - 0.077 J z 0. 157 J ± + 0.010 J z 0.348 0.124 

0.3 0.220(0.285) 0.293J± - 0.047 J z 0. 147 J± + 0.004 J z 0.297 0.101 

0.4 0.275 0.272J ± - 0.030J Z 0. 135 J± + 0.002 J z 0.253 0.083 

0.5 0.324(0.395) 0.250J± - 0.020J Z 0.124J ± + 0.001J Z 0.218 0.070 

0.6 0.370 0.230J± - 0.013J Z 0.114J X 0.188 0.060 
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FIG. 1. An example with a winding path showing the nonorthogonality of the approximate ground states |r) of the t — J z 
Hamiltonian. The configuration is generated by successively moving the hole in the perfect Neel configuration through the sites 
012301. But this is a configuration with a string of length 1 from site 2 to 1. 
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FIG. 2. The upper configuration has a hole separated from a pair of overturned spins by only one lattice site. The lower 
configuration has a hole separated from one isolated overturned spin by two lattice sites. Three hops of the hole in the lower 
configuration generate the upper configuration. 
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FIG. 3. This is a configuration with a string of length four starting from the site 1. Spin exchange between sites 1 and 2 
transforms it to a configuration with a string of length 2 starting from site 3. 
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FIG. 4. Same configuration as in Fig. |3j. Spin exchange between sites 2 and 3 transforms it to a configuration with a string 
of length 1 starting from site 4 and an isolated overturned spin at site 1. 



14 



